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Abstract 

We derive transport equations for fermions and bosons in spatially or temporally varying back- 
grounds with special symmetries, by use of the Schwinger-Keldysh formalism. In a noninteracting 
theory the coherence information is shown to be encoded in new singular shells for the 2-point 
function. Imposing this phase space structure to the interacting theory leads to a a self-consistent 
equation of motion for a physcial density matrix, including coherence and a well defined collision 
integral. The method is applied e.g. to demonstrate how an initially coherent out-of-equlibrium 
state approaches equlibrium through decoherence and thermalization. 
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1. Introduction 

Many problems in modern particle physics and cosmology require setting up transport 
equations for relativistic quantum systems in out-of-equilibrium conditions, including 
baryogenesis |ll2l3j , heavy ion collisions [1] and out-of equilibrium particle production [S] . 
It is straightforward to write a formal solution for the problem in the Shcwinger-Keldysh 
(SK) method, and detailed studies have been done on e.g. thermalization ^6, of quan- 
tum systems. It is much more difficult to find a simple enough an approximation scheme 
that is usable in practical applications. Here we report a recent progress |7|8j towards 
a scheme, in which one can account for the nonlocal coherence effects in the presence 
of hard collisions. The key observation is that in cases with particular symmetries (the 
homogenous time dependent case and stationary, planar symmetric case) the free 2-point 
Wightmann functions have new, hitherto unnoticed singular shell solutions (at fcp = 
for the homogenous system and at fc^ = in the static planar case), that carry the in- 
formation about the nonlocal quantum coherence. When this shell structure is fed into 
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Fig. 1. Dispersion relations in the case of a scalar mass function. The dark filled (red) blobs show the 
mass-shell contributions and the light (yellow) blobs show the corresponding coherence solution from 
the new fco = 0-shell in the homogenous case (left) and from the kz = 0-shell in the planar symmetric 
case (right). 

the dynamical equations for the full interacting theory, they reduce to a closed set of 
quantum transport equations for the weight functions of the singular shells. These func- 
tions correspond to the on-shell particle numbers (mass shells) and to a measure of the 
coherence (the new fco,z = 0-shells). We have demonstrated the correctness of our in- 
terpretation of the coherence shells through application to the Klein problem jE] ■ As an 
example of a case with collisions we have studied e.g. coherent production of unstable 
particles 8 . Here we will show how an initially highly coherent, out-of-equilibrium quan- 
tum field configuration approaches equilibrium through decoherence and thermalization. 
We will mainly be concerned with relativistic fermions, but the method is applicable also 
to scalar fields and nonrelativistic problems. 

2. Free theory and the shell structure 

Let us consider the SK-equations for a coUisionless fcrmion field, interacting with the 
background through a nontrivial mass term: 

C = lip ^-tjj + ■4'Lm{t, x)ipR + ipRm* {t, x)'iJjL ■ (1) 

In order to study the phase space structure of the correlator iG^{x,y) = {'ip{y)ip{x)) , we 
write the SK-equations in the mixed representation. In the coUisionless case, and to the 
lowest order in gradients the full SK-equation |3l8j reduces to {j(:+^^^—mii—imij^)G'^ — 
0. Consider first the homogenous case. Observing that helicity is a good quantum number, 
we can use the decomposition iG'f^j'^ = g'f^ ^{1 + hk ■ a) in the chiral representation. 
Equation of motion for then breaks into hermitean and antihermitean parts: 

(H): 2kog< = {H,g<}, (AH): idtg< = [H, g<] . (2) 

where H — —h\k\fj^~\-miiai — mia2 is the local Hamiltonian. The H-equation is algebraic, 
and it constraints the form of the possible solutions. Using Bloch representation 2g^ = 
g^ + g^ ■ (T, it can be written as Batg/} = 0. This equation has nontrivial solutions if 
det(i3) = (fcg — fc^ — |mp)A:Q = 0. Thus, somewhat surprisingly, in addition to the usual 
mass-shell solutions, there are new solutions at ko — (see Fig. ([T])). Explicit solutions 
on these shells are easily constructed: 



5h,m-s = 27rsfc(, f!}^^ (ko - h\k\(T3 + niRai - mia2)S{k^ - |m^|) 
g<„_, = {7r/\k\){ha3{mnf'^ - mjf^^) + |k|(ai/^ + a^f^)) 6{ko) , (3) 

where Skg = sgn(fco) and the mass- and coherence shell factors /^^ and are some 
unknown functions of t and k. A similar construction can be done for the stationary, but 
planar symmetric case |7]. In that case the spin s orthogonal to the symmetry plane is 
conserved, and eventually, in a special frame, one finds equations very similar to Eqs. Q: 

- 2sk,g< - Pg< + gf pt ^ ,,d^g< = pg< _ g< pt ^ (4) 

where P = kga^ + i{mfia2 + mjai) is the local mean field momentum operator. The 
algebraic equation again gives rise to a singular shell structure; only this time the mass- 
shells are accompanied by a singular shell at fc^ = [7j. The crucial issue leading to 
these new solutions is that we did not assume translational invariance in our treatment. 
If one imposes translational invariance, the dynamical equations in ^ and Q become 
algebraic as well and exclude the fcg^ = 0-solutions. 

The dynamical equations in Q and Q look compelling, but actually are meaningless 
because of the singular structure imposed by the algebraic constraints. Thus 5^j's are 
just phase space densities. Physical density matrices related to observable quantities 
can be introduced by convolving g^^'s with /cq- and/or fc^-dependent weight functions 
that quantify the prior extrenous information on the system [7]. Different types of such 
functions have been formally depicted in the right panel of Fig. [T] by dashed lines. In a 
typical application, say, to the early universe there is no prior information at all, and 
physically meaningful equations are obtained by integrating the singular equations of 
motion with a flat weight. 

The scattering off a step potential (the Klein problem) provides a nice test bench for 
our formalism because there, in the collisionless limit, equations ^ are exact equations 
for the problem. We applied eqs. Q to the Klein problem in ref. [7 , and our results 
demonstrated that the kz = 0-solutions are physical and describe the coherence: with 
them correct results for the transmission and reflection coefhcients and tunneling factors 
follows, while neglecting them completely destroys the quantum nature of the problem. 



3. Interacting fields 



With interactions included the full Schwinger-Keldysh equations become much more 
complicated of course. A tractable approximation scheme can still be obtained by com- 
bining the familiar quasiparticle approximation with our mean field limit treatment of the 
constraint equations. Eventually the approximation is that the constraint equations are 
corrected by inclusion of the real parts of the self energies, while the collision terms affect 
only the evolution equations. Note that these are just the standard assumptions that one 
makes when deriving the usual Boltzmann equation from the SK-equations j9]. Here, for 
simplicity, we neglect also the real part of the self energy. Furthermore, if we compute 
collision terms using thermal distribution functions, the generalization of eqs. ([2]) to the 
interacting case can be written in the following form: 

(H): 2kog<^{H,g<} 
(AH) : dtg< - ~t[H, g<] - {D, g< - {g<U} . (5) 
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Fig. 2. Evolution of the particle number = /+(|k|) (solid line) and the coherence fc = (/^ + /|)^^'^ 
(dashed line) for an initially highly coherent out-of-equilibrium state under decohering interactions [S]. 

where the precise form of the coUision term D depends on the interactions. When the 
evolution equation is integrated over the fco-variable with a flat weight, contributions 
from each of the shells depicted in the left panel of the Fig. 2 contribute to the physical 
density matrix p. The evolution equation can be written in the form 

dtPh = -iiH, ph] - Ig{ph) , (6) 

where the commutator term describes the usual quantum mechanical evolution of the 
density matrix, and the integral Ig describes the effect of collisions. (For an explicit form 
of Ig for a particular example see ref. [8].) The four independent components of ph are 
in one-to-one correspondence with the on shell functions fj^, a = ±,1,2. This allows 
for the collision integral to be evaluated in terms of ph , leading to closure as is formally 
indicated in eq. ([6| . It is important to note that the collision integral Ig is nontrivial; it is 
computed at off shell momentum for the coherence shell. As a sample application of our 
method, we show in fig. [2] how an initially highly coherent out-of-equilibrium quantum 
state relaxes to equilibrium as a result of collisions [8] . 

Our method can also be applied to the bosonic fields with similar results. The coher- 
ence again resides at fco,z = 0-shells and a closed set of moment equations that include 
coherence and collisions can be found [10 . It is also easy to write down similar quantum 
transport equations in the case of nonrelativistic fields |10j . 
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